Guided Vortex Motion in Superconductors with a Square Antidot Array 
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We have measured the in-plane anisotropy of the vortex mobility in a thin Pb film with a square 
array of antidots. The Lorentz force, acting on the vortices, was rotated by adding two perpendicular 
currents and keeping the amplitude of the net current constant. One set of voltage probes was used 
to detect the vortex motion. We show that the pinning landscape provided by the square antidot 
lattice influences the vortex motion in two different ways. First, the modulus of the vortex velocity 
becomes angular dependent with a lower mobility along the diagonals of the pinning array. Second, 
the vortex displacement is preferentially parallel to the principal axes of the underlying pinning 
lattice, giving rise to a misalignment between the vortex velocity and the applied Lorentz force. 
We show that this anisotropic vortex motion is temperature dependent and progressively fades out 
when approaching the normal state. 



The interaction of a moving elastic medium with an 
array of obstacles has received much attention during the 
last years, in part due to a variety of condensed matter 
systems that can be described within this model. 1-4 In 
general, the symmetry of the substrate over which the 
elastic medium moves, should also be reflected in the 
dynamic properties of the system. 

Particular interest has gone to the study of the flux 
line lattice (FLL) in superconductors with an artificial 
pinning potential, such as randomly distributed colum- 
nar tracks, 5 arrays of submicron holes (" antidots 1 '' ) 6-8 or 
magnetic dots. 9 ' 10 

In a superconductor with a random distribution of pin- 
ning centers, an isotropic vortex mobility is expected 
with a vortex motion parallel to the Lorentz force. On 
the other hand, if the rotational symmetry of the pinning 
potential is broken, the response will be anisotropic, 11 
and guided vortex motion may appear. 

Recent experiments 12 have shown that the presence 
of twin boundaries in high temperature superconductors 
leads to a preferential motion of the FLL along the twin- 
ning planes. Guided vortex motion has also been found 
in systems with a four-fold symmetry, like square arrays 
of Josephson junctions. 13 Theoretical work on vortex mo- 
tion in a square array of pinning centers has forecasted 
guided vortex motion or channeling along the principal 
axes of the array, and along its diagonals. 14 " 17 In these 
systems, the vortex velocity v preferentially snaps to the 
symmetry axes of the pinning array, leading to a mis- 
alignment between vortex velocity and the applied driv- 
ing force f . Even though this lock-in transition has been 
predicted in many numerical simulations for a square ar- 
ray of pinning centers, it had not yet been observed ex- 
perimentally. 

In this work, we report on electro-transport measure- 
ments on a Pb thin film with a square array of antidots, 
for the whole 360° range of in-plane current orientations. 



This allows us to determine the misalignment between 
the applied current and the electric field generated by the 
vortex movement. We show clear evidence of guided vor- 
tex motion along the principal axes of the antidot array, 
and determine the trapping angle in the studied struc- 
ture. 

The used sample is a 50-nm-thick Pb film, with a 
square antidot lattice (antidot size b 2 = 0.6 x 0.6/zm 2 ) of 
period d = 1.5fj,m, which corresponds to the first match- 
ing field Hi = &o/d 2 =9.2 G, with $ the superconduct- 
ing flux quantum. We patterned the Pb film in a cross- 
shaped geometry (see Fig. 1) to allow electrical transport 
measurements at any in-plane orientation of the applied 
current. The cross consists of two 300 /im wide strips 
containing the square array of antidots aligned with the 
principal axes of the cross. 
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FIG. 1. Left: sketch of the cross-shaped Pb film with a 
square antidot array (grey area). Right: zoom in of the central 
part of the cross and orientation of the total applied current 
Jt and vortex velocity v with respect to the x-axis of the 
pinning array. The black dots represent the voltage probes. 

This pattern was prepared by electron-beam lithog- 
raphy in a polymethyl metacrylate/methyl metacrylate 
(PMMA/MMA) resist bilayer covering the Si0 2 sub- 
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strate. A Ge(20 A)/Pb(500 A)/Ge(200 A) film was then 
electron-beam evaporated onto this mask while keeping 
the substrate at 77 K. Finally, the resist was removed in 
a liftoff procedure in warm acetone. The sample has a 
T c =7.220 K with a transition width of 7 mK for H=0. 

The transport measurements were performed in a 3 He 
cryostat using two independent dc currents, Jx and 
Jy, applied at the extremes of the cross legs, allow- 
ing to control the orientation (<fi) of the total current 
Jt = Jx cos0 x + Jy sin</> y, with respect to the prin- 
cipal axes of the square pinning array. The voltage con- 
tacts were aligned with the x direction (see Fig. 1, right 
panel). In the Vx(4>) measurements, we sweep the an- 
gle cf) of the total c urrent Jt while keeping its amplitude 
Jt = V 1 J\-V Jy constant. The magnetic field H was 
applied perpendicular to the film surface. 

It is worth to notice that in previous works, 12,18 in 
order to determine the direction of the vortex move- 
ment, two sets of voltage probes arranged perpendicu- 
larly (aligned with each of the current directions), have 
been used. The simultaneous acquisition of the indepen- 
dent components of the electrical held E = v x H, allows 
one to determine the modulus and direction of the av- 
erage vortex velocity v and therefore the misalignment 
between v and the applied Lorentz force f = Jt x 3*0- 
Although this procedure is strictly necessary for uni- 
axial systems, simple symmetry considerations associ- 
ated with the periodic pinning landscape show that for 
higher degree of symmetry, it is possible to obtain the 
same information with a single pair of voltage contacts. 
Indeed, in our experiment, the intrinsic four-fold sym- 
metry imposed by the pinning potential implies that 
Vy(4>) = Vx(4> + 90°), and therefore by measuring Vx(4>) 
in the whole angular range, we can deduce Vy (</>). 

In Fig. 2(a) we show the field dependence of the re- 
sistance R(H) for J T = 3kA/cm 2 and T =7.20 K, for 
several angles <f> of the total current. The most obvious 
features of this figure are the well defined dips at \/2Hi 
and Hi, as a consequence of the reduction of the FLL 
mobility, reminiscence of the vortex configurations ener- 
getically stable in the equilibrium state. For low vortex 
density (H < OAHi), a linear field dependence charac- 
teristic of a flux flow regime ppp — p n H/H C 2, with p n 
the normal state resistivity and H C 2 the upper critical 
field, is observed. 19 

It should be noted that, since the experiments were 
carried out very close to T c , temperature stability plays 
a crucial role. For instance, the curve at 4> = 50° in 
Fig. 2(a) exhibits a small asymmetry between the positive 
and negative branches of the magnetoresistance (R(H) ^ 
R(—H)) because of a temperature drift of only 0.2mK. 




FIG. 2. Magnetoresistance R(H) at Jt= 3kA/cm 2 for (a) 
several angles cf> of the total current Jt at T = 7.20 K (b) 
4> = 0° and T = 7.17 K. In (b), the field dependence of the 
eccentricity factor e = V((j> = 45°)/V(0 = 0°) is shown (open 
symbols). 

Due to the misalignment between the wire bounded 
voltage probes and the lithographically defined current 
contacts, the resistance in the normal state reaches the 
condition R = (i.e. Vx = 0) at <fi = 77° rather than at 
(f> = 90°. Remarkably, a finite resistance appears when 
penetrating in the superconducting phase, \H\ < 15 G 
« l.QHi (open cicrcles in Fig. 2(a)). Since the R(H) 
curve at </> = 77° is symmetric around H = 0, we rule 
out the possibility to attribute this dissipation to the Hall 
effect. We have also checked that thermoelectric voltages, 
generated by thermal gradients, have a negligible contri- 
bution to the total voltage. This "even Hall effect" might 
be attributed to guided vortex motion, 12,18 however we 
have observed that this effect persists in plain films with- 
out any periodic structure, where transport properties 
should be isotropic. Similar anomalies have also been 
reported in the transverse voltage in Josephson Junction 
arrays 20 and in the longitudinal voltage in plain films. 21 
Although there is no general consensus about the origin 
of this effect, a possible explanation involves proximity 
effects due to the N/S boundary in the voltage contacts. 22 

Fig. 2(b) shows R{H) for = 0° and J T = 3fc Aj cm 2 
at T =7.17 K. For fields \H\ < Hi, where vortices are 
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strongly pinned at the antidots, the dissipation lies be- 
low our experimental resolution. At this temperature, 
the maximum number n s of flux quanta that an anti- 
dot can hold is n s w 6/4£(T) - 0.3 < l, 23 where £(T) 
is the superconducting coherence length. Under these 
circumstances, interstitial vortices appear in the sample 
for magnetic fields H > H\. They are weakly pinned 
( "caged" ) at the interstitial positions between the repul- 
sive saturated antidots. The matching features at 
and 2Hi are a result of their stable vortex configurations. 
In this scenario, where two species of vortices coexist, the 
motion of interstitial vortices should be strongly affected 
by the presence of those vortices which remain pinned. 
As field progressively increases, vortex-vortex interaction 
becomes more relevant in comparison with the vortex- 
pinning energy and eventually, close to the upper critical 
field H c2l a vortex liquid phase appears. In order to 
study the role of the periodic pinning array on the vor- 
tex dynamics in each of these regimes, we have measured 
the in-plane angular dependence of the voltage Vx{4>) 
for the same conditions of Fig. 2(b) {Jt — 3kA/cm 2 and 
T =7.17 K) at H =10, 12, 14, 15, 16, 18, 20 and 22 G. 
The curves Vx {4 1 ) for some of them are presented in Fig. 
3. 
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FIG. 3. Angular dependence of the voltage Vx{4>) for 
Jt = ZkA/cm 2 and T =7.17 K, at H =10, 18 and 22 G. 

The upper panel of Fig. 3 shows Vx{4>) close to the 
normal state (H c2 (7.17K) w 25 G). The recorded signal 
follows a smooth sine function dependency characteristic 
of an isotropic regime with a constant vortex velocity in- 
dependent of the current orientation. At H —18 G, near 
to the second matching field, the response is quite dif- 
ferent. Now, a sawtooth like curve clearly manifests the 



existence of an in-plane anisotropy in the vortex move- 
ment. This anisotropy becomes even more evident as H 
approaches to the first matching condition (lower panel) . 
As we will show below, the characteristic shape of Vx {<fi), 
alternating well pronounced peaks and plateaus, is a fin- 
gerprint of guided vortex motion. 

As we pointed out previously, knowing the longitudi- 
nal component of the voltage Vx (</>), we can deduce the 
transversal component Vy(4>) = Vx(</>+90°). This allows 
us to build a polar graph Vx x Vy , as shown in the inset of 
Fig. 4 for H =10, 14, 16 and 22 G, where the curves have 
been scaled for the sake of clarity. The circular shape of 
the outer curve (H =22 G) reflects the rotational invari- 
ance of an isotropic conductor. As field decreases, the 
shape of the polar plot smoothly evolves from a circle, 
through a rhombus (not shown) and eventually data lies 
along an astroid for \H\ < 20 G, indicating that vortices 
move faster along the principal axes [±1,0], [0, ±1] (easy 
direction) of the pinning structure than along the diago- 
nals (hard direction). 
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FIG. 4. Angle 6 of the vortex velocity as a function of the 
Lorentz force orientation a for Jt = 3fc A/ cm 2 and T =7.17 
K, at H =10 G. Inset: polar graph Vx xVy for Jt = 3kA/cm 2 
and T =7.17 K, at (from inner to outer) H =10, 14, 16 and 
22 G. 

The degree of anisotropy of the vortex velocity is given 
by the eccentricity factor e = V(<f> = 45°)/F(0 = 0°). In 
case of a circular shape, e = 1, whereas for a rhombus 
s = l/ In Fig. 2(b) we show the anisotropy e from the 
analysis of the data at T=7.17 K, as a function of field. 
We clearly observe that e{H) reproduces the behavior of 
the magnetoresistance R(H). According to the theory, 13 
the vortex dynamics should become more anisotropic as 
Jt decreases towards the depinning current Jd p - On the 
other hand, at constant Jt, a higher Jd p yields a lower 
dissipation. Therefore, a stronger anisotropy should be 
detected where; R(H) has local minima, in agreement 
with our observations. A similar e{H) behavior, but 
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somewhat less pronounced was found at higher temper- 
atures for all fields \H\ < H c2 . This result indicates that 
the anisotropic vortex motion is present regardless of the 
dynamic regime. 

Once the polar graph is built, it is possible to extract 
the angle of the electric field 9e = Atan{yy jVx) and 
hence the angle of the vortex velocity 9 = 6e + 90° as a 
function of the Lorentz force orientation a — <j> — 90°. In 
the main panel of Fig. 4 we show 9(a) so obtained for 
H =10 G, folded to the fourth and first quadrant. For 
an isotropic conductor, 9 = a and a straight line should 
be observed. In contrast, a highly non-linear response is 
found. 

Starting from a = and rotating the current in 
counter-clockwise direction, i.e. increasing the angle <j>, 
we observe that for low values of a, the force f lies ahead 
of v, which remains close to the [1, 0] axis (9 ~ 0°). The 
misalignment \a — 6\ reaches its maximum value at the 
lock-in angle ~ 15°, beyond which the vortex move- 
ment is no longer restricted to follow the [1,0] direction 
and a fast increase in 9 is observed. At the [1, 1] orienta- 
tion, force and velocity are aligned due to the four-fold 
symmetry of the underlying pinning potential. Above 
[1,1], the Lorentz force lags behind the velocity vector 
which rapidly approaches to the [0, 1] axis. This is the 
central result of the paper which clearly indicates the role 
played by the periodic pinning potential in the vortex dy- 
namics. 

Recent simulations 14 for a moving FLL in a supercon- 
ducting films with square pinning array show that minor 
lock-in phases should also be present at ±45° orienta- 
tions. Further features might be present near the driving 
angles along the symmetry directions of the pinning array 
other than the principal axes. 15 However, no evidence of 
such phases were found in our experiments. The absence 
of these regimes may be attributed to the relatively large 
size of the antidots. Indeed, since the pinning potential 
is averaged in the direction of motion, it is expected that 
the movement at ±45° (diagonals of the square antidots) 
be progressively washed out as the size of the antidots in- 
creases. Well defined channels will exist only if d/Ab > 1. 
In our experimental conditions d/4b = 0.625 and thus 
channeling effects should be suppressed. This assump- 
tion is confirmed in the weakly coupled superconducting 
wire network limit by simulations for square Josephson 
Junction arrays by Marconi et al. where no particular 
effects appear at ±45°. 

In summary, using dc-transport measurements we have 
demonstrated that the interaction of a moving flux line 
lattice with a periodic pinning potential, breaks the ro- 
tational symmetry and induces a guided vortex motion 
along the principal symmetry orientations of the pinning 
array. As a consequence, we found that for orientations 
different than the symmetry directions, the vortex ve- 
locity is not aligned with the Lorentz force. This effect 
vanishes approaching to the normal state as the vortex 
movement becomes insensitive to the pinning symmetry. 

We would like to thank Clecio C. de Souza Silva for 



helpful discussions. This work was supported by the 
ESF "Vortex" Program, the Belgian Interuniversity At- 
traction Poles (IUAP), the Flemish GO A and FWO Pro- 
grams. 



1 P. Martinoli, Phys. Rev. B 17, 1175 (1978). 

2 K. Runge and B. Pannetier, Europhys. Lett. 24, 737 
(1993). 

3 A. van Blaanderen, R. Ruel, and P. Wiltzius, Nature 385, 
321 (1997). 

4 J. Hu and R. M. Westervelt, Phys. Rev. B 55, 771 (1997). 

5 D. Niebieskikwiat, A. Silhanek, L. Civale, G. Nieva, P. 
Levy, and L. Krusin-Elbaum, Phys. Rev. B 63, 144504 

(2001) . 

6 M. Baert, V. V. Metlushko, R. Jonckheere, V. V. 
Moshchalkov, and Y. Bruynseraede, Phys. Rev. Lett. 74, 
3269 (1995). 

7 V. V. Moshchalkov, M. Baert, V. V. Metlushko, E. Rosseel, 
M. J. Van Bael, K. Temst, R. Jonckheere, and Y. Bruynser- 
aede, Phys. Rev. B 54, 7385 (1996). 

8 V. V. Moshchalkov, M. Baert, V. V. Metlushko, E. Rosseel, 
M. J. Van Bael, K. Temst, Y. Bruynseraede and R. Jonck- 
heere, Phys. Rev. B 57, 3615 (1998). 

9 D. J. Morgan and J. B. Ketterson, Phys. Rev. Lett. 80, 
3614 (1998). 

10 M. J. Van Bael, K. Temst, V. V. Moshchalkov and Y. 
Bruynseraede, Phys. Rev. B 59, 14674 (1999). 

11 L. Van Look, B. Y. Zhu, R. Jonckheere, B. R. Zhao, Z. 
X. Zhao, and V. V. Moshchalkov, Phys Rev. B 66, 214511 

(2002) . 

12 H. Pastoriza, S. Candia, and G. Nieva, Phys. Rev. Lett. 
83, 1026 (1999). 

13 V. I. Marconi, S. Candia, P. Balenzuela, H. Pastoriza, 
D. Dommguez, and P. Martinoli, Phys. Rev. B 62, 4096 

(2000) . 

14 C. C. de Souza Silva and Gilson Carneiro, Phys. Rev. B 
66, 054514 (2002). 

15 C. Reichhardt and F. Nori, Phys. Rev. Lett. 82, 414 (1999). 

16 Gilson Carneiro, Phys. Rev. B 62, R14661 (2000); Phys. 
Rev. B 66, 054523 (2002). 

17 V. I. Marconi and D. Dominguez, Phys. Rev. B 63, 174509 

(2001) . 

18 Vicent Berseth, Ph.D. Thesis N°2017, Ecole Politechnique 
Federale de Lousanne, Switzerland (1999). 

19 E. Rosseel, M. Van Bael, M. Baert, R. Jonckheere, V. 
V. Moshchalkov, and Y. Bruynseraede, Phys. Rev. B 53, 
R2983 (1996). 

20 B. J. van Wees, H. S. J. van der Zant, and J. E. Mooij, 
Phys. Rev. B 35, R7291 (1987). 

21 Y. K. Kwong, K. Lin, M. S. Isaacson, and J. M. Parpia, 
Phys. Rev. Lett. 65, 2905 (1990); M. Park, M. S. Isaacson, 
and J. M. Parpia, Phys. Rev. Lett. 75, 3740 (1995). 

22 M. Park, M. S. Isaacson, and J. M. Parpia, Phys. Rev. B 
55, 9067 (1997). 

23 GS. Mkrtchyan and V.V. Schmidt, Zh. Eksp. Teor. Fiz. 61, 



4 



367 (1971) [Sov. Phys. JETP 34, 195 (1972)]; H. Nordborg 
and V. M. Vinokur, Phys. Rev. B 62, 12408 (2000). 



5 



